Tc-Enhanced Codoping Method for GaAs-based Dilute Magnetic Semiconductors 



m 
o 
o 

(N 
> 

o 



C/2 



o 



> 



o 



I 

O 

o 



Van An DINH^-^Q Kazunori Sato^, and Hiroshi Katayama-Yoshida^-^ 
^Department of Condensed Matter Physics, and ^Department of Computational Nanomaterials Design, 
Nanoscience and Nanotechnology Center, The Institute of Scientifie and Industrial Research, 
Osaka University, 8-1 Mihogaoka , Ibaraki, Osaka 567-0047, Japan. 
^ Institut fiir Festkorperforchung, Forschungszentrum Jiilich, D-52425 Jiilich, Germany 

Based on ab initio calculations of Gai-ajMn^^NyAsi-y and Gski-xMiixCyAsi-y, we propose a new 
codoping method to enhance the Curie temperature Tc of diluted magnetic semiconductors. The 
solubility of Mn can be increased up to high concentration by the codoping of N or C to reduce 
the lattice and volume expansion caused by Mn doping. It is found that the impurity band of the 
majority spin is strongly broadened and pushed up into the higher energy region due to the strong 
p-d hybridization caused by the codoping, and the Tc becomes higher than the room temperature 
at x> 6%. 
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The discovery of ferromagnetic diluted magnetic 
semiconductors 111 H (DMS) by incorporating the 3d 
transition metal atom (TM) into the III-V compound 
semiconductors was a great step towards the application 
of the spin degree of freedom of the electrons in semi- 
conductor spintronics. One of the most important issue 
regarding the realization of the semiconductor spintronics 
is to synthesize the DMS with a Curie temperature (Tc) 
higher than the room temperature. The III-V compound- 
semiconductor-based DMS such as (Ga,TM)As is one 
of the most promising candidates from the view-point 
of the industrial application. Recently, the ab initio 
calculations^ ^3 have predicted the high Tc in Mn or 
Cr-doped GaAs. However, up to now, the Tc which was 
reported in the new experiment[3 of (Ga,Mn)As is limit- 
ted to 150K because of the low solubility of TM atoms 
in GaAs. 

The Tc of the ferromagnetic (Ga,Mn)As is expected to 
increase with both the Mn and hole concentrations P, 0, 
H^ . In fact, the Mn atom at the substitutional site of 
the GaAs acts as an acceptor, but in (Ga,Mn)As, the 
hole concentration is substantitially lower than the Mn 
concentration. It is ascribable that when the TM atoms 
such as Mn, Cr, etc. are doped into GaAs, the lattice 
constant and volume of the (Ga,TM)As become larger 
than that of the host material, leading to the presence of 
compensating donors formed by interstitial TM atoms. 
As pointed out by Yu et al. Q, the Mn atom can oc- 
cupy the tetrahedral interstitial sites in the zinc-blende 
structure and acts as a double donor. The compensation 
of the substitutional Mn acceptors by the interstitial Mn 
donors decreases the hole concentrations and Tc. We can 
reduce the acceptor compensations due to the interstitial 
Mn donors by reducing the lattice or volume expansion 
using the codoping method. Nitrogen (N) or Carbon (C) 
is one of the excellent candidates for this purpose by us- 
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ing the codoping of Mn and N (or C) since the N and C 
have the smaller atom radius than the As atom. 

Since N 2p-level is lower than the As 4p-lcvel, we can 
expect strong p-d hybridization in Mn-N than Mn-As. 
Therefore, the antibonding Mn 3(i- impurity state (ta) is 
pushed up into the higher energy region in the band gap. 
Based on the double-exchanged mechanism of the ferro- 
magnetism in DMS, Tc is propotional to ^/x and n^d{Ep) 
[T3II14 [l5l |. where x is the concentration of 3c?-TM, and 
nzdiEp) is the 3d-TM partial density of states (PDOS) 
at the Fermi level Ep- It is well known that the band 
structure of GaAs changes strongly with the incorpora- 
tion of isoelectronic N at As sites to form GaNAsp. At 
N concentration of 10%, the GaNAs is expected to be- 
come semi-metal0. Therefore, besides the avoidance of 
the substitutional Mn acceptor compensation by the in- 
terstitial Mn donors, an influence of strong p-d hybridiza- 
tion caused by N codoping on the magnetic properties of 
(Ga,TM)(N,As) to increase the n-^diEp) and Tc is also 
expected. 

Following the above-mentioned idea of materials de- 
sign, we propose a new codoping method to enhance the 
Tc of GaAs-based DMS, by which TM and N (or C) are 
codoped simultanously into Ga and As sites of GaAs, re- 
spectively. Incorporated concentrations of N (or C) in 
the codoping method are chosen to satisfy the condition 
of no lattice expansion or no volume expansion. 

When GaAs is doped by Mn atoms, the lattice con- 
stant a of the host material changes, and the dependence 
of the deviation of the lattice constant on the Mn concen- 
tration obeys the Vegard law[3|. Moreover, the volume 
of the Mn-doped GaAs becomes also larger. On the con- 
trary, when some As atoms of the GaAs are substituted 
by N (or C) atoms, the lattice constant and volume de- 
crease. In order to reduce the Mn donors at tetrahedral 
interstitial sites and increase the solubility of substitu- 
tional Mn acceptors, we set the N concentration y so 
that the lattice constant a and volume V are kept un- 
changed based on the experimentally observed Vegard 
law. In other words, the dependence of j/ on 2; is chosen 
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so that the deviation of the lattice constant Aa{x,y) or 
the volume expansion of the lattice cells AV{x, y) can 
be ignorable. As reported by Fan et al.^^, the change 
of the lattice constant in the case of N doping with low 
concentration(?/ < 3%) still obeys the Vegard law. There- 
fore, we can assume here that the deviation of the lattice 
constant is linear with N concentration. The lattice con- 
stants of 5.6533 A, 4.52 A and 5.98 A for GaAs, GaN and 
MnAs, and the atom radii of 1.221 A, 1.367 A, 1.245 A 
and 0.549 A for Ga, Mn, As and N, respectively, are used 
to derive the dependences of the N fraction on the Mn 
concentration. We obtained the following relations for 
the options of the codoping method based upon the ex- 
periments U,E^:y = 0.29x (for Aa = 0), ?/ = 0.42x (for 
AV = 0). 

The calculations are based on the KKR-CPA-LDA 
method 0. The MACHIKANEYAMA-2000 package 
coded by H. Akai (Osaka University) is used for 
the calculation of the density of states (DOS) and the 
Tc of Gai_a;Mna;NyAsi_,y. In the present method, the 
Mn atoms replace randomly the Ga atoms in GaAs 
host material. The As atoms are substituted ran- 
domly by the N atoms. The Tc^s evaluated by us- 
ing the mean field approximationj^ ITsl [l3 | . For com- 
parison, the calculation of the of the (Ga,Mn)As 
is also performed. To simulate the spin-glass state of 
Gai_a;Mn2;NyAsi_j,, the Ga cations are substituted ran- 
domly by Mn^,j.f2 and Mn^^f2 magnetic ions, where 
t and J, denote the direction of the local moment 
of Mn atoms. Hence, the ferromagnetic and spin- 
glass states can be decribed as (Gai_2^,Mn1^j;)(N.y,Asi_.y) 
and (Gai_3;,MnT^/2,Mn''-^/2)(Nj^,Asi_j^), respectively. 
Throughout the calculations, the potential form is re- 
stricted to the niufhn-tin type, and the mufhn-tin radii 
of 1.224 A are chosen for both cations and anions. The 
wave functions in each muffin-tin sphere are expanded 
with real harmonic up to / = 2, where / is the angular 
momentum at each site. The relativistic effect is also 
taken into account by using scalar relativistic approxi- 
mation. 328 fc-points in the irreducible part of the first 
Brilouin zone arc used in our calculations. 

Figure H depicts the DOS of Gai_j;Mnj;NyAsi_y for 
several concentrations of Mn {x = 2%, 4%, 8%, 12%). 
The anti-bonding states ta of 3d-Mn caused by the hy- 
bridization with 4p-states of As ions and 2p-states of N 
ions form an impurity band in the band gap. With incre- 
scasing Mn concentration, the N concentration is also ac- 
cordingly increased upon the codoping, this anti-bonding 
impurity band ta connects with the host valence band to 
widen the band width (W) which is proportional to the 
root square of the Mn concentration {W oc ^/x) , and the 
band gap is considerably narrowed (see Fig. ^d)). Such 
a behavior relates to the well-known band gap bowing 
effect in the GaAs incorporated with N ^Jl- Moreover, 
the Fermi level of the majority spins shiftes into the new 
part of the valence band to realize a half-metallic fcrro- 
magnctism. 

It is well known that the stabilization of the fcrro- 
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FIG. 1: Total DOS per unit cell (solid curve), PDOS of 3d- 
states per Mn atom at Ga site (dashed curve) and 2p-state per 
N atom at As site (dotted curve) for several Mn concentra- 
tions. N concentration is set to be 29% of Mn concentration. 
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FIG. 2: fl-PDOS of Mn 3d-states in 

(Gai_a;,Mn^i)(Nj;,Asi_j;). N is incorporated at As site 
with concentration being 29% of Mn concentration. 



magnetism in GaAs-based DMS is due to the double- 
exchange mechanism 0, 0, . In the DMS with the 
double-exchange mechanism being dominant, the Tc is 
propotional to the band width [Tc (x W oa ^/x) and the 
PDOS of the anti-bonding states ta at the Fermi level 
(Tc (X ny,d{Ep)). In the present case, the ta-peaks are 
sensitive with the change of Mn concentration. The im- 
purity band width W is proportional to the Mn concen- 
tration: W cx which is shown in Fig. ^ and also in 
Fig. El 
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FIG. 3: PDOS of Mn 3d-states for x = 5%. Incorporated 
N concentration is chosen as 0% (dashed curve), 29% (sohd 
curve) and 42% (dotted curve) of Mn concentration. 
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The dependence of the ta-band width of majority spins 
on the Mn concentration varying from 3% to 15% is 
shown in Fig. El Similar to the Fig. Q the N concen- 
tration is chosen to be equal to 29% of Mn concentration 
to compensate the lattice expansion upon the codoping. 
It is easy to see that ta-bands are gradually broadened 
with the substitutional Mn concentration. Since the sta- 
bilization energy of the fcrromagnctism is proportional 
to the to-band width, it suggests to open up a realistic 
possibility to synthezise a high- Tc DMS based on GaAs 
host material by the codoping of Mn and N. 

For comparison, we ploted also the PDOS 
of the Bd-Mn states of the ferromagnetic state 
(Gai_a;,Mn1^^)(Nj^,Asi_j^) at 5% of Mn concentration of 
the both codoping and single-doping cases(Fig. The 
dashed curve refers to the 3(i-PD0S of Mn ions in single- 
doped GaMnAs (i.e. y ~ Q). The remainders illustrate 
the Sd-PDOS for two options of the N codoping method. 
It is remarkable that at 5% of Mn concentration, the 
3(i-PDOSs of Mn atoms in the codoping cases show 
a large difference from that in the single-doping case. 
Because of the host valence band originates from anion 
p-states, the doping of N with 2p-level deeper than the 
shallow 4p-level of As into the host GaAs leads to the 
change of the energy levels of the anti-bonding ta-statcs 
with strong p-d hybridization, and ta-states become 
more 3(i-character. By the codoping of Mn and N, the 
peaks of Sd-Mn PDOS of fo-states are pushed up and 
closer to the Fermi level than that in the single-doping 
case. The ia-peaks are pushed into the higher energy 
region in the band gap, and the ferromagnetism caused 
by the double exchange interaction is more stablized due 
to the higher density of 3d-Mn PDOS at Ep- 

It is emphasized here that the codoped N plays an 
important role to increase the Tc in Gai^j^Mn^jNyAsi-j,: 



FIG. 4: Curie temperature vs. Mn concentration of 
Gai-ajMuiiAs (crosshair-dashed curve), Gai-aiMniNyAsi-y 
and Gai_a;Mna;Cy Asi-y (solid curve), y varies as 0.291 
(triangle-solid curve), 0.42a; (circle-solid curve) for the codop- 
ing of N, and 0.33a: for the codoping of C. The dotted line 
indicates the room temperature (300 K) 



(i) Smaller atomic radius of N atoms compensates the 
lattice expansion caused by Mn doping at the Ga sites, 
then increases the solubility of Mn atoms in GaAs. (ii) It 
changes the band structure of the host material, resulting 
in the shift of the 3d- states of ta towards the higher 
energy region in the band gap, the increase of the PDOS 
of ta-states at the Fermi level, and therefore the increase 
of the Tc caused by the enhancement of the ferromagnetic 
double exchange interaction. 

As usual, the estimation of the Tc is performed by us- 
ing a mapping on the Hciscnberg model in the mean 
field approximation ^1_3. 14]: fc^Tc = 2Ai?/3x, where 
/S.E is the total energy difference between the spin-glass 
and ferromagnetic states. Fig. 0] shows the comparison 
of the Tc between the single-doped Gai_a;Mna;As, and 
codoped Gai_2;Mni:Nj,Asi_y. As seen from Fig. 0] the 
Tc increases with subtitutional concentration as \/x due 
to the ferromagnetic double-exchange interaction. The 
Tc of the single-doped GaMnAs gains the room tem- 
peratures at 10% of Mn concentration and reaches the 
saturation regime at a; w 12%. However, the Tc of 
Gai_a;Mn2;Nj,Asi_j, gains the room temperature at the 
smaller Mn concentrations {x > 5% for the doping op- 
tion y = 0.42a; and x > 6% for y = 0.29a;), and has 
a tendency to increase (~ 30%) with the whole of Mn 
concentrations. At the higher Mn concentrations, the 
codoping method gradually enhances the Tc due to the 
presence of N atoms. 



4 



It is noted that while the highest Tc of Gai-ajMn^^As 
was reported about 110 K by Matsukura et aZ.0] and 
newly 150 K by Ku et al.^, the ab initio calculations 14, 
Ha . [iJI predicted the much higher Tc . This discrepancy 
is caused by the compensation effect in the experiment 
since the thermal equilibrium solubility of Mn atoms in 
GaAs is very small 10^^ cm~'^). This effect can be re- 
duced by the codoping method as proposed in this letter, 
by which we can reduce the interstitial Mn atoms by N 
codoping with reducing the lattice expansion caused by 
Mn doping. However, it is noted here that since the influ- 
ence of the band gap bowing effect due to N doping, the 
band gap becomes narrower, resulting in the ia-states 
could be located in the valence band, then the p-d ex- 
change interaction becomes dominant at the sufficiently 
high concentration of the N atoms. Therefore, the codop- 
ing with N concentration higher than 10% may causse the 
reduction of the Tc- 

Since the 2/3 of the ij^-bands are occupied by electron 
in (Gai_a;,Mnl^a;)(Nj,,Asi_y) as shown in Fig. 2, we can 
further increase the by the p-type codoping which in- 
creases n^diEp) and the stabilization energy from the 
ferromagnetic double-exchange interaction. Carbon is 
suitable candidate for this purpose. Using the ion radii 
of 0.62 A, 0.67 A, 0.58 A, and 0.16 A for Ga, Mn, As 
and C, respectively, we derived the relation of Mn and C 
concentrations as y = 0.33a; regarding no volume expan- 
sion. Our ab initio calculations show the more gradual 
enhancement of the Tc compared with the codoping of N 
as demonstrated in Fig. 4. 

In conclusion, we propose a new method of codoping to 
synthezise the High- Tc DMS based on the ah initio calcu- 



lations of the electronic structure and Curie temperatures 
of the Gai_2;Mna;NyAsi_y and Gai-^MnajCyAsi-j,. Mn 
and N (or C) are incorporated into host materials on 
condition that they do not change the lattice constant or 
volume of the host materials to increase the solubility of 
substitutional Mn atoms, as well as, the PDOS of the t\ 
states at the Fermi level. 

The obtained results for Gai_a;Mna;NyAsi_j, and 
Gai_2;Mn2:CyAsi_y show the gradual enhancement of 
the Tc by the codoping method. It is suggested that 
since we can reduce the compensating donors formed by 
the interstitial Mn atoms by reducing the lattice expan- 
sions upon the codoping, Tc could be higher than the 
room temperature at a; > 6% for the codoping of N, and 
X > 5% for the codoping of C. The influence of the N 
atoms on the magnetic properties with deeper 2p— levels 
of N and strong p-d hybridization between Mn 3d- and 
N 2p- levels is also strongly expected to increase the Tc- 
Also, the p-type codoping of C is the promising method 
to gradualy enhance the Tc based on the ferromagnetic 
double-exchange interactions. 
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